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Abstract 

For a particular set of Boltzmann weights and a particular boundary con- 
dition for the six vertex model in statistical mechanics, we compute explicitly 
the partition function and show it to be equal to a factorial Schur function, 
giving a new proof of a theorem of Lascoux. 

1 Introduction 

A state of the six vertex model, also called square ice, is defined as follows. 
Between each pair of neighbouring vertices of (a subset of) the square lattice 
in two spatial dimensions, we place an oriented edge, subject to the condi- 
tion that at each vertex, exactly two of the four edges are oriented inbound. 
Around each vertex there are ( 2 ) = 6 potential configurations, hence the name. 
We shall always consider states of the six vertex model that are bounded spa- 
tially with fixed boundary conditions. One can also realise this model in terms 
of a square array of ice molecules, where an oxygen atom is placed at each 
vertex, and a hydrogen atom is placed on each edge, with the covalent bonding 
determined by the hydrogen atom being bonded to the source oxygen atom 
of its orientation. It is the ice interperetation that we shall use in this paper. 
In diagrams, we represent oxygen atoms by circles and hydrogen atoms by 
smaller solid disks, with the existence of covalent bonds being exhibited via 
spacial proximity in diagrams. 

To each vertex v in a state T, we associate a Boltzmann weight w v (T) 
which depends only on the location of v in the lattice and the local configu- 
ration of incoming and outgoing edges containing the vertex v. The partition 
function (from statistical mechanics) for this model is then defined to be 

T v 

where the sum is taken over all possible states T and the product is over all 
vertices v in the model. It is this function which we intend to evaluate in this 
paper. 

Let n be a positive integer and A a partition with n parts. So A = 
(Ai, A2, . . . , A„) where Ai > • • • > \ n > are all integers. We consider a 
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rectangular array of square ice consisting of n rows and n + Ai columns (the 
effect of adding more columns to the right will trivially have no effect on 
the partition function). On this model, we impose the following boundary 
conditions, which are the same as found in [BBFJ. The left and right hand 
boundaries will be filled with hydrogen atoms with the lower boundary com- 
pletely devoid of hydrogen atoms. These conditions already imply that along 
the upper boundary exactly n of the columns will not have hydrogen atoms, 
we insist that these be in the columns Aj + n + 1 — i for i = 1, . . . , n. We shall 
call this boundary condition the A-boundary condition. An example of such 
a configuration with A = (5, 4, 1) 



is 



Let xi, X2, ■ ■ ■ , x n and oj, 02, . . . , a n+ \ 1 be two alphabets of variables. For 
computing the partition function, we use the following set of Boltzmann 
weights. For the vertical molecule in row i and column j, we attach the 
weight Xi/aj. For the molecule with hydrogen atoms to the north and west, 
we attach the weight Xj/a,- — 1. For all other orientations of molecules, we 
attach the weight 1. Let Z\(x\a) denote the corresponding partition function. 
Our main theorem is the following. 



Theorem 1.1. \Las01/\ Theorem 1] Up to a monomial, the partition function 
Z\ is equal to a factorial Schur function. Explicitly, we have 



Z\(x\a) = -p—yS X (x\a). 



The notation of the above, used through out this paper, is that whenever 
y = (yi, y2, ■ ■ •) and fi = (pi, H2, • • •) is a multi-index, then y^ = y^y^ 2 
We define the partition p to be (n — 1, n — 2, . . . , 0) and the multi-index 5 as 
(0, 1, . . . , n — 1). For any partition A, A' denotes its conjugate (the reflection of 
A, identified with its Young diagram, in the main diagonal). A cell € A 
refers to the box in the intersection of the i-th row and j-th column of the 
Young diagram of A. 

Our proof is fundamentally different to that of [Las07j. We begin by 
showing the symmetry of x~ s Z\(x\a) and then use the characteristic vanishing 
properties of factorial Schur functions to identify the partition function as a 
factorial Schur function. 

The author would like to thank Ben Brubaker, Alain Lascoux and Steven 
Sam for helpful conversations. 



2 Factorial Schur functions 



We present an overview of the definition and basic properties of factorial Schur 
functions for which |Mol98] can be taken as a reference. 
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Definition 2.1. The factorial Schur function s\(x\a) is defined to be 

S\(x\a) = ^ Y\.( X T(a) ~ a T (a)+c(a)) 
T aGA 

where c(a) is the content of the cell a in the Young diagram of \, defined by 
c(i,j) = j — i. The summation is taken over all semistandard Young tableaux 
T of shape A. 

We now state the primary properties of these factorial Schur functions 
without proof. 

Proposition 2.2. The factorial Schur function s\(x\a) is symmetric in the 
variables x%, . . . ,x n . 

Proposition 2.3. The collection of factorial Schur functions s\(x\a) as A 
runs over all partitions with n parts, forms a basis for the ring of symmetric 
polynomials in x\, . . . ,s n . Restricting to symmetric polynomials of degree at 
most m, the collection of factorial Schur functions s\(x\a) as above with the 
additional restriction that |A| = Y2i ^« — m f orms a basis for this vector space. 

Given a partition p = . . . , /%), define the sequence a M = ((a M )i, . . . , (a M ) n ) 
by («fy)t = a n+ i_ i+ A r 

Proposition 2.4 (Vanishing Theorem). Suppose that A and p are partitions 
of at most n parts. Then 

s\{a^\a) =0 if \(£ /j,, 

s x (a\\a) = Y[ (a-n+i-i+x, - a n _ A / +i ). 
(M')eA 



3 Strict Gelfand-Tsetlin Patterns and Stair- 
cases 

To understand the possible states of square ice with the A-boundary condition, 
and to connect our combinatorial setup with the staircase language of Lascoux 
[Las07| . it will be useful to consider a bijection between such states and two 
other combinatorial objects, namely strict Gelfand-Tsetlin patterns with top 
row A + p and staircases with rightmost column missing A + p. We now define 
these other families of objects. 

A Gelfand-Tsetlin pattern is a triangular array of natural numbers 

tl,l *1,2 *1,3 • • • H,n-1 H,n 

*2,1 h,2 ^2,3 • • • *2,n-l 

% = \ . . > 

subject to the inequalities tij < U+xj < tjj+i for all i and j. This pattern is 
further called strict if the additional inequalities tij < ijj+i also hold. 

To say that the top row of X is equal to A + p is to say that ti jn+ i_,; = 
n + 1 — i + Aj for all i. 
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A staircase is a semistandard Young tableau of shape (Ai + n, X± + n — 
l,...,Ai)' filled with numbers from {1,2, ...,Ai + n} with the additional 
condition that the diagonals are weakly decreasing in the south-east direction 
when written in the French notation. 

To say that a staircase has rightmost column missing A + p means that 
the integers in the rightmost column comprise all integers from 1 to Ai + n 
inclusive with the exception of those integers of the form n + 1 — i + Aj for 
some i. 

An example of a staircase for A = (5, 4, 1) is 
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We have chosen to present this particular example because under the bijection 
we are about to construct, it gets mapped to the sample array of square ice 
given in the introduction. 

Proposition 3.1. Let A be a partition. There are natural bijections between 
the following three sets of combinatorial objects 

1. States of the six vertex model with the A boundary condition, 

2. Strict Gelfand- Tsetlin patterns with top row A + p, and 

3. Staircases whose rightmost column is missing A + p. 

Proof. For each state of the six vertex model, and each % = 2, . . . , n, there are 
n + 1 — i occurrences of an oxygen atom between rows i — 1 and i covalently 
bonded to the hydrogen atom in row i — 1 above it. If we let tij for j = 
1, . . . , n + 1 — i be the column numbers of these oxygen atoms, then this 
corresponding collection of integers forms a strict Gelfand-Tsetlin pattern 
with top row A + p. 

Given a strict Gelfand-Tsetlin pattern T with top row A + p, we construct 
a staircase whose rightmost column is missing A + p in the following manner: 
We fill column j + 1 with integers u%, . . . , Un-j^^ such that 

{l,2,...,n + Ai} = {m,... } U {tn+l-j,l, ■ ■ 

It is easily checked that these maps give the desired bijections. □ 

With these bijections in hand, it is a straightforward manner to check that 
our formulation of Theorem 11.11 is equivalent to |Las071 Theorem 1] , noting 
that a double Schubert polynomial for a Grassmannian permutation is the 
same as a factorial Schur function. 
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4 The star triangle relation 



We consider a set of Boltzmann weights on a diagonally oriented vertex, cor- 
responding to string numbers i and j where i runs from lower left to top 
right. 

For each of the six admissible orientations of the ice molecule, we attach 
the following weights: 



1 • 1 • • f 1 1 

- , # Q : - , # Q # : 

— s X% "^3 — \ ^3 

Now consider the following diagrams 



O :0 - 



o 



1 

Xj_ 



1 



e<5 





We consider boundary conditions a, (3, 7, S, e and £, exactly three of 
which correspond to the placement of a hydrogen atom in the corresponding 
place in the above diagrams. With these boundary conditions, we consider 
all possible ways of completing this to an arrangement of square ice by the 
placement of exactly one hydrogen atom along each string marked with a 
question mark. Each such state has a corresponding Boltzmann weight, where 
for the diagonal molecules we attach the weights presented immediately above 
while for the rectilinear molecules, we attach the same Boltzmann weight as 
in the introduction, where the row number is determined by the string it is 
sitting on. 

Theorem 4.1 (Star Triangle Identity). The star triangle identity holds in the 
following sense. We fix boundary conditions and sum the Boltzmann weights 
of all possible states of the left hand diagram. Then this equals the sum of all 
the Boltzmann weights of all possible states of the right hand diagram with the 
same boundary conditions. 

Proof. There are a total of ( 3 ) = 20 possible boundary conditions for which 
this theorem needs to be checked. We will illustrate this with a particular 
example. 

Consider the boundary condition consisting of a hydrogen atom in the /?, 
7 and £ positions. Then our purported identity becomes 




Xi 1 



(---)- + -(--iH = H-i) 

X i X j CL X j Ob Ob CL (X X j 

The other 19 cases (the majority of which are simpler than this) are left to 
the reader. □ 
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Corollary 4.2. We have the identity 
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Z\{. . . ,Xi,Xi+i, ...\a) = —Z\(. . . ,x i+ i,Xi, . . . 



a). 



Proof. The left hand side is the partition function of our square ice model, 
with a crossing between horizonal strings i and i + 1 added on the left. All 
boundary conditions are the same as for the calculation of Z\. The right 
hand side is the partition function for the same ice model, but this time 
the crossing between strings i and % + 1 is added on the right. Using the 
star triangle relation proved above, we can incrementally push the crossing 
of horizonal strings past each vertical string while preserving the partition 
function, hence proving the Corollary. □ 

Corollary 4.3. The function x~ s Z\(x\a) is symmetric in x%, . . . ,x n . 

Corollary 4.4. The function x~ 5 Z\{x\a) is a polynomial in x%, . . . ,x n . 

Proof. A priori, we know this to be a Laurent polynomial in x\,...,x n that 
is a polynomial in x\. Since it is symmetric, it must thus be a polynomial in 

x\ , . . . , x n . n 

We remark that the symmetry result proved in this section does not use the 
upper or lower boundary conditions imposed on our states of square ice, and 
thus gives a direct proof of the symmetry property of the function F{u, v; z) 
from |Las07} Theorem 2]. 



Our method of proof is to show that the function Z\{x\a) possesses the same 
characteristic vanishing properties of the factorial Schur functions. 

Theorem 5.1. For two partitions A and // of at most n parts, we have 



Proof. Fix a state, and assume that this state gives a non-zero contribution to 
the partition function Z\(a^\a). Under the bijection between states of square 
ice and strict Gelfand-Tsetlin patterns, let ki be the rightmost entry in the 
i-th row of the corresponding Gelfand-Tsetlin pattern. We shall prove by 
descending induction on i the inequality 



For any j such that fcj+i < j < k{, there is a factor (xi/aj — 1) in the Boltzmann 
weight of this state. We have the inequality n + 1 — i + ^ > n + 1 — (i 1 + 1) + 
fJ-i+i > ki + i by our inductive hypothesis. Since (a^)i = a n+ i_j + Ai, in order 
for this state to give a non-zero contribution to Z\{a^\a), we must have that 
n + 1 — i + fii > ki, as required. 



5 Proof of the Main Theorem 



Z\(a^\a) = unless AC/i, 




n + 1 - i + [ii > ki. 
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Note that for all i, we have ki > n + 1 — i + \ . Hence \ii > \ for all i, 
showing that /i D A as required, proving the first part of the theorem. 

To compute Z\(a\\a), notice that the above argument shows that there 
is only one state which gives a non-zero contribution to the sum. Under the 
bijection with Gelfand-Tsetlin patterns, this is the state with ij = t\j for all 
The formula for Z\(a\\a) is now immediate. □ 

Now we have collated enough data to be able to prove Theorem II .li 

Proof. We can see, for example using the bijection with Gelfand-Tsetlin pat- 
terns, that Z\ is a polynomial of degree at most [A + p| in x. Thus Z\/x s is 
a symmetric polynomial of degree at most |A| in x, so we may write it in the 
form. 

x~ s Z\ = c^(a)s^{x\a) 

£t:|/^|<|A 

for some coefficients which are rational functions of the a variables only. 

Inducting on fx, we substitute x = a M into the above equation to conclude 
that c M (a) = for all \i ^ A with < |A[. By substituting x = a\, we are 
also able to calculate c\(a) = (a^ x+ ^ completing the proof. □ 
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